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Abstract. Several algorithms are available in the literature for finding the 
entire set of Pareto-optimal solutions in MultiObjective Linear Programming 
(MOLP). However, it has not been proposed so far an interior point algo- 
rithm that finds all Pareto-optimal solutions of MOLP. We present an explicit 
construction, based on a transformation of any MOLP into a finite sequence 
of SemiDefinite Programs (SDP), the solutions of which give the entire set 
of Pareto-optimal extreme points solutions of MOLP. These SDP problems 
are solved by interior point methods; thus our approach provides a pseudo- 
polynomial interior point methodology to find the set of Pareto-optimal solu- 
tions of MOLP. 

1. Introduction 

Although aheady more than 60 years old, Linear Programming is still nowa- 
days one of the most important areas of research and application in Mathematical 
Programming/Operations Research. There are many different algorithms for its 
resolution but essentially all of them can be classified in four main branches: pri- 
mal simplex, dual simplex, primal-dual simplex and interior point methods. 

Although not as the same level as standard linear programming, it is also com- 
monly accepted that Multiobjective Linear programming (MOLP) is another very 
important area of activity within the optimization field. Multiobjective optimiza- 
tion is motivated by the need to consider multiple, conflicting, objectives in real 
world decision making problems. Multiple objective linear programming has been 
a subject of research since the 1960s for its relevance in practice, as a mathematical 
topic in its own right. Its development has come in parallel to the scalar counterpart 
and its theory and algorithms have been developed mainly in the last two decades 
with some remarkable exceptions. Solving a MOLP means to obtain the entire set 
of Pareto-optimal solutions. By analogy with the scalar case, one can find in the 
specialized literature several algorithms to find the entire Pareto-optimal set. In 
the multiobjective case there are primal simplex- like type algorithms (Steuer |26j . 
Yu and Zeleny |29j and the references therein) , primal-dual simplex- like algorithms 
(Ehrgott, Puerto and Rodriguez-Chia [5]) and dual simplex-like algorithms (Benson 
[2], Ehrgott, Lohne and Shao [S]). Moreover, there are some partial approaches that 
use interior point methods to approximate or to find some Pareto-optimal points 
(see e.g. Fliege [UlIE]) but no interior point algorithm that finds all Pareto-optimal 
solutions has been proposed so far. Quoting [F': "Interior point algorithms either 
require interaction with a decision maker to find a preferred efficient solution or 
find at most an efficient face. No interior point algorithm that finds all efficient 
solutions has been proposed" . 

Due to the natural parallelism between these two areas, namely scalar and mul- 
tiobjective linear programming, different authors wondered whether there would 
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exist a fully based interior point approach to generate the complete Pareto-optimal 
set of a MOLP. This paper gives an affirmative answer to this question. 

The goal of this paper is to develop an interior point method that generates the 
entire Pareto-optimal set of MOLP. The interest of this research is two-fold: 1) 
It strengthens the parallelism between these two areas of Mathematical Program- 
ming, namely Linear Programming and Multiobjective Linear Programming; 2) It 
is theoretically appealing because shows how to adapt some tools available nowa- 
days in the field of polynomial optimization ([Ml [17]) to be applied in a completely 
different field as is Multiobjective Linear Programming. 

From the results in this paper, we show that it is theoretically possible to find the 
entire set of Pareto-optimal solutions using interior point algorithms. Indeed, we 
will explicitly describe an algorithm to perform that task. Our approach consists of 
constructing a finite sequence (its cardinality m is the number of constraints in the 
MOLP) Semidefinite programs based on the original MOLP, the solution of which 
can be used to generate the entire Pareto-optimal set using the moment matrix 
algorithm by Lasserre |16| . Nevertheless, this program size is not polynomial in 
the input size of our original MOLP. There is no surprise in this fact because it 
is well-known that finding the entire Pareto-optimal set is NP-hard since it may 
be equivalent to enumerate the vertices of its feasible region Khachygan et al [5]. 
Thus if our construction were polynomial in the input size, interior point algorithms 
would prove polynomiality of the problem. Needless to say that even though our 
construction is explicit, we do not claim that this approach is computationally 
competitive with other approaches available nowadays, as for instance the variant 
of the outer algorithm of Benson developed recently in [8]. Nevertheless, it is of 
theoretical interest because proposes the first interior point based algorithm for 
obtaining the entire Pareto-optimal set, it suggests a completely different scheme 
to address this problem and proposes another pseudo-polynomial algorithm to deal 
with that computation. It is worth mentioning a different framework where a 
similar approach is applicable: multiobjective polynomial integer programming. 
We point out that the use of algebraic tools for solving multiobjective is not new. 
The interested reader is referred to ^j4, for further details. 

The rest of the paper is organized as follows. Section[2|describes the general mul- 
tiobjective optimization problem, together with the considered concept of solution 
and it recalls the main results that are needed for the rest of the paper. Section [3] 
presents the theoretical results that lead to obtain the entire set of Pareto-optimal 
solutions of MOLP using interior point algorithms. Here, we explicitly describe a 
system of polynomial equations that encodes the entire set of Pareto-optimal ex- 
treme points which is the basis of our next results. Section [4| reduces the problem 
of finding the entire set of Pareto-optimal extreme points to solving m explicit SDP 
problems. They are feasibility problems since their objective functions are constant. 
In addition, we also show how to extract all its, finitely many, feasible solutions by 
applying the moment matrix algorithm 16J. This construction is illustrated with 
an example taken from the literature [8 . In the final section of the paper we draw 
some conclusions. 

2. Multiobjective linear programming 

In this section we recall the main theoretical results for the development in this 
paper. We begin by describing the general framework to cast the problem to be 
handled. A multiobjective optimization problem (MOP) consists of: 



(MOP) 



?;-min {fi{x), . . . , fk{x)) 

s.t. x e 5 
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where fi : M" — > M for i — 1, . . . ,/c are the objective functions and S C M" is 
the feasible region. The symbol v — min means that we want to minimize aU the 
objective functions simultaneously. If there is no conflict between the objective 
functions, then, a solution can be found where every objective function attains its 
optimum. In such a case, no special methods are needed. Otherwise, first we need 
to state what we understand by a solution of the above problem. It is commonly 
agreed that the solution concept is related with the notion of Pareto-optimal points. 
The definition of Pareto-optimal solution is due to Edgeworth [7]. However, that 
name was first used by Koopmans in [TF after the developments by Pareto in 



Definition 1. A decision vector x* ^ S is a Pareto-optimal solution for (MOP) 



if it does not exist another decision vector x S such that fi{x) < fi{x*) for all 
i = 1, . . . , k and fj{x) < fj{x*) for at least one index j . 

The set of Pareto-optimal solutions is called the Pareto-optimal set. (PO set, 
for short.) 



If X* is a Pareto-optimal solution f{x*) is said an efhcient solution of (MOP). 



It is commonly agreed that solving a (MOP) consists of finding the entire set of 
Pareto-optimal solutions. There are some other solution concepts for multiobjective 
optimization problems as local, weak, proper or strong Pareto optimality (see |20j). 
In this paper we will restrict ourselves to the standard definition of Pareto-optimal 
solutions, although similar approaches, to the one adopted in this paper, would be 
also valid for the rest of the solution concepts in multiobjective optimization. 

One of the main methods for describing the Pareto-optimal set of a multiobjective 
optimization problem is by scalarization, that is, transforming the multiobjective 
problem into a single or a family of single-objective problems with a real- valued ob- 
jective function, depending on some parameters. This enables the use of the theory 
and methods of scalar optimization to be applicable to get the solutions of MOP. 
The importance of these methods rests on the fact that Pareto-optimal solutions 
of MOP can be characterized, in most cases, as solutions of certain single objective 
optimization problems. There are several of these scalarization methods for solving 
MOP (see [2D])- Among them, we consider the weighting method. The idea of this 
method is to associate each objective function with a weighting coefficient so as to 
minimize the weighted sum of the objective functions. The weighting method can 
be used so that the decision maker specifies a weighting vector representing his pref- 
erence information. However, this method can also be used to generate iteratively 



several solutions of (MOP ) by modifying adequately the weighting coefficients. The 
success of this approach is based on the following result by Gass and Saaty [13] or 
Zadeh [5D]. 

Lemma 2. Let fi be convex for all i = 1, . . . ,fc and S be a convex set. Then, 



if X* £ S is a Pareto-optimal solution of (MOP), there exists a weighting vector 
A e M^\{0}, Ai = 1 such that X* is a solution of the following scalar problem: 

k 

(SP) min^A,/,(a;) 

1=1 
s.t. X £ S 

Note that from the above result, in particular if fi is linear for all i = 1, . . . , fc 



and 5 is a convex polyhedron, all the Pareto optimal solutions of (MOP) can be 
found by the weighting method. 

In this paper we are interested in solving a special class of multiobjective op- 
timization problems where all the objective functions are linear and the feasible 
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region is described by a set of linear inequalities, that is: 
(MOLP) V - min Cx (c^x, c'^x) 

s.t. Ax>b 
x>0 

with the i-the row of C € i = l,...,k,Ae M"^" and b E M™. We assume 

w.l.o.g. that the system Ax > b, x > 0, has not redundant inequalities and that 
^ cone{c^ , . . . , c'^) since otherwise all the feasible region is Pareto-optimal and 
the problem is trivial. Under our assumption, the Pareto-optimal set is included in 
the boundary of the feasible region and it is well-known that it is edge connected 
but not necessarily convex, see e.g. [26]. With this settings we say that a Pareto- 



optimal solution for ( MOLP I is an extreme point if it is a vertex of the polyhedron 



defining the feasible region of ( MOLP ) 



In addition and in order to simplify our presentation, we will assume w.l.o.g. that 
the feasible region is a polytope and that we are given redundant upper bounds on 
the values of the x variables, namely we are given ubj" such that xj < ubj', i = 
1, . . . ,n. Note that by the fact that the feasible region is a polytope these bounds 
can always be obtained for sufficiently large ub^ values and they are redundant. 

Lemma [2] ensures that to solve (MOLP I it suffices to apply the above weighting 
method. In doing that, this problem is transformed to a family of parametric linear 
programming problems. For this reason, we recall here some results about linear 
programming that will be useful in the next sections. 

Consider the following pair of dual linear programming problems: 



(LP) 



mm c X 
s.t. Ax>b 
x>0 



max 6* 



(DLP) 



s.t. A^u < c 

u>0 



with A e 



be 



and c e 



By our assumption on the feasible region of the primal problem ( LP I , we observe 



that the dual problem can be also assumed to be bounded and that we can assume 
that we are also given redundant upper bounds on the feasible values of (DLP). 
Namely, we know ubf such that Ui < ubf for all i = 1, . . . , m. 

The following classical result, whose proof can be found in [55], gives the rela- 
tionship between the optimal solutions of the above problems: 



Lemma 3 (Strong Duality Theorem). Let x* be a feasible solution of (LP) and 
let u* be a feasible solution of (DLP I such that cx* = b^u* , then x* is an optimal 



solution for (LP) and u* is an optimal solution for (DLP I 



An important consequence of the Strong Duality Theorem is the following result, 
usually called Complementary Slackness Theorem. It is useful in order to detect 



whether a feasible solution of ( LP I and a feasible solution of ( DLP ) are optimal 



of their respective problems. Its proof, that can be found in i24j, is easily deduced 
from Lemma m 

Lemma 4 (Complementary Slackness Property). Let x* be a feasible solution of 
(LP I and let u* be a feasible solution of (DLP). Then, the following statements 



are equivalent: 



(1) X* is an optimal solution of (LP) and u* is an optimal solution of (DLP|. 

(2) X* and u* satisfy u*\b - Ax) = and (u*A - c^)x* = 0. 

In the next section we will show how the entire set of Pareto-optimal solutions of a 



multiobjective linear problem ( MOLP ) can be obtained by solving m-SDP problems 
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being m the number of non-redundant inequalities of the original problem. This 
shows that interior point methodologies can be also used to determine the Pareto- 



optimal set of (MOLP) 



3. A POLYNOMAL SYSTEM OF INEQUALITIES ENCODING THE PARETO-OPTIMAL 

SET OF MOLP 



Note that by Lemma [2] solving ( MOLP ) is equivalent to solve the following 
parametric family of single objective problems: 

fe 

(LP a) min^A, 



\ic'''x 



1=1 
s.t. Ax> h 



for all A G A := {A G 



!}■ 



For each A e A, the dual of (LPj]) is: 



(DLPa) 



max Uj bj 



s.t. u*A < ^ Aic' 



■u > 

Hence, by Lemma |4j a solution of (MOLP) must be a solution of the following 
system of polynomial equations/inequalities: 

u\b-Ax) = 
fc 



(Sysl) 



{J2 - u*A)x = 
Ax>b 



i^A < ^ A, 

1=1 



u,x > 0, A e A 

Solutions of this system are triplets (x, u, A) such that x is a Pareto-optimal solution 



of (MOLP), and optimal solution of {PL\); and u is optimal solution of {DLP\). 



We will call such a triplet a 'valid triplet'. However, we observe that the above 
system may have an infinite number of solutions because there may be a continuum 
number of solutions in x and u since the Pareto-optimal set is a connected union 
of faces of the polyhedron Ax > b, x > 0. Nevertheless by the edge-connectedness 
of the Pareto-optinal set, it suffices to know the Pareto-optimal points that are 
extreme points in the feasible region to reconstruct the entire set (see for instance 
[1112]). Therefore, our goal is to reduce the solutions of the above system to a finite 
number of the original ones (extreme points) that are enough to reconstruct the 



entire set of Pareto-optimal solutions. Hence, we shall transform System (Sysl ) so 
as to characterize only extreme points of the original feasible region. 

It is worth mentioning that a similar approach could be applied to general mul- 
tiojective convex optimization. Indeed, in that case we have to replace the use 



of System (Sysl) by using Kuhn- Tucker optimality conditions on the paramet- 
ric family of weighted problems. The only inconvenient is that in that case the 
system may have a continuum of solutions and in general, one cannot build the 
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entire Pareto-optimal set with a finite set of representatives, as we do in ttie linear 
case. Tliis implies that finiteness results cannot be obtained. Another different 
framework where a similar approach is applicable is in multiobjective polynomial 
combinatorial optimization. The interested reader is referred to [31 0] for further 
details. 

As it is usual, let B denote a basis of A, i.e. a full rank submatrix of A. In 
addition and when there is no possible confusion, we will also use B as the set of 
the indices of its columns. Analogously, N denotes the set of columns of A not 
in B and are the coefficients of the ^-th objective function that corresponds 
to variables in the basis B. Finally, with xb and xn, we refer to the variables 
corresponding to the columns in B and N , respectively. 

Let us denote by (U-Sys-B-il for B a full rank submatrix of A and i — 1, . . . ,m, 
the system: 

k 

(U-Sys-B-i) 



Xic^ - u*_^A^, > A,. 



E 



Xic'sB-^A,j - ^ic] < 0,Vj G N 



5: A, 



1 



where U-i, A^i stand, respectively, for the vector of variables u without Ui and the 
matrix A where we have removed the i-th row. In the same way, Ai. is the z-th 
row, A.j is the j-th column and Aij is the (i, j) element of A, respectively. 

Finally, let M be the least common multiple of all the determinants of full rank 
submatrices of A and Mi be the least common multiple of all the determinants of 
full rank submatrices of (U-Sys-B-i), for all z = 1, . . . , to and B. 



Theorem 5. If x is a Pareto-optimal solution, extreme point of the feasible region 



of (MOLF) then Mx is the projection onto the first n-components of a solution of 
at least one of the systems (Sys — il, for some i — 1, . . . ,to; 

fe 

(1) h°{\) :=^A,-M, = 



hi{x, u) 



u\Mh- 

k 



Ax) = 



h{{x,u) X(C^ - ut,A_i - MiA,.)x = 

£=1 

(2) g°{x) := As.x - Mb, > 0, s = 1, . . . ,to, 

fe 

(3) g]{u, A) := ^ A,c^^ - uU{A^i); - M,A,, > 0, j = 1, . . . , n, 

£=1 
ub^ M 

Pj{x):= J|(a;j-£)=0, j = l,...,n, 
i;(A):-[](A,-£)=0, r = l. 



1=0 
ub°Mi 

(us-^) = 0, s = 1, . . . ,m, s 7^ I, 



1=0 
Mi 



(Sys - i) 



i=o 
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Conversely, any of the finitely many solutions of the systems (Sys — i| for all 
i = l,...,m induces a Pareto- optimal solution of (MOLP) and all the Pareto- 
optimal extreme points are included among them. 



Proof. Let x be a Pareto-optimal extreme point solution of (MOLP). Being an 



extreme point of the feasible region Ax > h, x > Q means that there exists a 

basis B which defines x such that x = ^ "'(^ ^ """" ^ n-i? 

MxB 


by Lemma [2j x must be optimal for the problem {PLx[x)) for some A(a;) (Note 
that this weighting coefhcient may depend on x). The optimality condition of x 
translates into the following necessary and sufficient condition for any valid \{x): 



clearly Mx — 



and XB = B^^b > 0. Next, 
, has integer coordinates in the range [O^ub^M]. Besides, 



(4) 



k 

E 



Xic'bB^^A.. 



k 



V j e N, 



(5) 



where as usual N denotes the set of columns of A not in B and Cg are the coefficients 
of the ^-th objective function that correspond to variables in the basis B. 

Moreover, by Lemma[4] each pair {x, \{x)), where X{x) is defined by Q-ljsj), has 
associated an optimal extreme solution of {DLPx(x)) such that {x , u{x) , X{x)) is a 
valid triplet, i.e. it satisfies System (Sysl). By the fact that u{x) is an extreme 
solution of {DLPx(x)) it must satisfy: 



(6) 



u*A 



k 

E 

1=1 



\ic^ < 



Assume w.l.o.g. that for some valid triplets {x , u{x) , X{x)) the i-th component of 
u{x) is positive, i.e. u{x)i > 0, then there is always one of these valid triplets for 
X so that {u{x),X{x)) is an extreme solution of the system of linear inequalities 
defined by Q, ([5| and This system written in matrix form is; 



X 



e 

AU 




CbB 



c 
1 



c. 



N 



X 



< 
< 





At 

1 



Next, any extreme solution of this system must be associated with a basis of the 
matrix Li. Therefore if Mi is the least common multiple of the determinants of all 
full rank submatrices of Li, then (Mi(w_i, 1), Af,;A) has integer coordinates in the 
range [0,m5^M,] for u and [0,Mi\ for A. 

Finally, it is clear by our construction that {Mx, Miu{x), MiX{x)) is a solution 

of[S^O 

Conversely, by Lemma |4j any solution (x , u{x) , X{x)) of Sys — i for some i = 
l,...,m defines a Pareto-optimal triplet {x/M,u{x)/M„X{x)/Mi) of (|MOLP[). In 



addition, we have proved above that all Pareto-optimal extreme point solutions of 

for i 



( MOLP I are among the solutions of the systems ( Sys — i 



concludes the proof. 



1, 



, m which 
□ 



The above transformation makes use of big upper bounds to ensure that extreme 
points of all systems of rational inequalities that come from feasibility and optimal- 
ity of valid triplets are integer. In most cases, those bounds can be strengthened 
taking advantage of the particular structure of the problems given rise to sharper 
bounds. In particular if the primal, the dual or both are integer polytopes one can 
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take M = 1, M, =^ 1, for alH = 1, . . . , m or both M = = 1 for alH = 1, . . . , m; 
respectively. Eventually, we will need only to transform the range of the lambda 
variables to make them integer. 

To finish this section, we briefly remark the meaning of Theorem [6j We have 
proved that the entire set of Pareto-optimal extreme points is encoded in the set of 
solutions of TO system of polynomial inequalities, projecting the first n components 
of their solutions. (Observe that this set may include some extra Pareto-optimal 
solutions not being extreme points.) The finiteness of the solution sets of these 
systems is crucial to develop an exact method to get the entire Pareto-optimal set, 
as we will show in the next section. 



4. Semidefinite programming versus Multiobjetive Linear 

programming 

In this section we describe how to obtain the Pareto-optimal set of any MOLP 
by applying tools borrowed from the theory of moments and SDP. We use standard 
notation in the field (see e.g. [HI). 

We denote by IR[a;, u, A] the ring of real polynomials in the variables x = (xi, . . . , a;„), u — 
(ui, . . . , Um), A = (Ai, . . . , Afe); and by Mfx, u, X]d C M[a;, u, A] the space of poly- 
nomials of degree at most d G N (here N denotes the set of nonnegative inte- 
gers). We also denote by ;B = {^"u'^A''' : (a,/3,7) £ l^"x™x'^} a canonical basis of 
monomials for R[x, u, A] , where x°'u^\^ — a;"^ • • • x'^"u^^ ■ ■ ■ A^^ • • • A^*" , for any 
(a,/3,7) e N"^"^^ 

For any sequence indexed in the canonical monomial basis B, y — (yap-y) C M, let 
Ly : M.[x, u, A] — M be the linear functional defined, for any / = X]Q07gN"X"x'» /a/57 x°'u^X^ G 

R[a;,U, A], as Ly(/) := X]a/37eN"X™x'' /a/37 ya;37- 

The moment matrix Md(y) of order d associated with y, has its rows and columns 
indexed by {x'^u^X'') and Mrf(y)(a/37, a'/3'7') := Ly(a:("+"')u(^+^') A(''+''')) = 

?/(„+c.')(/3-H/3')(7+7')' l"'/3'7'l < d. 

For g := XiieNixmxfe g&{xu\Y € K[a;,M, A], the localizing matrix yij_{g^y) of or- 
der d associated with y and has its rows and columns indexed by {[xuXY) and 
Mrf(5,y)(a/?7,a'/3V) Ly (a;("+"')u(^+^') A(^+^')5(a;wA)) = 552/5+(a+a')(/3+/5')(7+7')' 

for |a^7|, |a'/3'7'| < d. 

Depending on their parity, let 2C,j or — 1 be the degree oipj, j = 1, . . . ,n; 2r]\ 
or 2r7j — 1 be the degree of , s = 1, . . . , to, i = 1, . . . , to, s ^ i and or 2vl — 1 
be the degree of , r — 1, . . . , fc, i — 1, . . . , to. Recall that these functions were 



defined associated with some of the constraints that appear in ( Sys — i I . Finally, 
for any symmetric matrix P by P ^ we refer to P being positive semidefinite. 
With this notation, we are in position to present our next result: 



Theorem 6. The entire set of Pareto-optimal extreme point solutions of (MOLP) 
can be obtained by solving m semidefinite programs {SDP'' — N* ) i = 1, . . . ,m, for 
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some N* e N. 
{SDP' - N*) 



aim 2/0 ■— 1 
s.t. Mn.{x,u,\) >z 0, 





MAr._i(/lO(A)) = 0, 






MN'-i{h'i{x,u))=0, 






MN.^i{hl{x,u,X))^0, 








,m, 




MN'-i{9}{u,Xj) hO, j^l,. 


. . ,n 


(7) 


Mn._^^{Pj{x))^0, j = 


. ,71, 


(8) 


M^._^.(g:(u)) =0, 5 = 1,.. 


. ,n. 


(9) 


Mw._,.(t;(w)) = 0, r = l,.. 


.,k. 



Proof. The orem [5] proves that the entire set of Pareto-optimal extreme point so- 
lutions of (MOLPJ) can be obtained as the projection of the solutions of systems 
( Sys — i I, i = 1, . . . , m. Each of these systems is defined by a closed, compact semi- 
algebraic set with finitely many feasible solutions that by compacity, clearly satisfy 
Putinar's condition (see [16]). Hence, we can apply [16, Theorem 6.1] to conclude 
that there exists N* < +oo so that all the sol utions of Sys tem ( Sys — i ) can be ob- 

m. □ 



tained by solving the feasibility SDP problem ( SDP^ — N 



for alH = 1 , 



We note in passing that the finiteness of N* is explicit and it can be bounded 
above by some known constants which depend on the input data. The interested 
reader is referred to [23] and [5S]. 

As a direct consequence of the equivalence between MOLP and solving a finite 
series of m-SDP (Theorem [6| ; and the fact that SDP can be solved by interior 
points algorithms we conclude the following corollary. 



Corollary 7. All Pareto-optimal extreme solutions of (MOLP) can be obtained 
using interior point algorithms. 

It is well-known that solving SDP problems is polynomially doable. Therefore, 
one may conclude from the above reformulation that obtaining the entire Pareto- 



optimal set of (MOLP I is polynomial. However, we cannot conclude this result 
from Theorem |6 The reason being that the dimension of the m-SDP problems that 
need to be solved are not polynomial in the size of the original problem because N* 
is finite but its size may be exponential in n, to, k. Therefore, the above approach 
only gives a pseudo-polynomial approach to obtain the Pareto-optimal set of MOLP 
programs. 

Finally, we want to describe, based on the above results, an explicit methodology 



to enumerate the real solutions of the above closed semi-algebraic sets ( Sys 
i = 1, . . . , m. (Note that this is equivalent to enumerate all the optimal solutions of 
a polynomial optimization problem where the objective is a constant; see Theorem 



In order to do that, first we transform w.l.o.g. (Sys — i) into an algebraic set, 
in the standard manner, by simply adding non-positive slack variables w € R™, 
z G K" to the constraints ^ and ([S]), respectively. Thus, from now on we assume 

are in equation form. 



that all the constraints in the SDP relaxations ( SDP 

For any i = 1, . . . , to, consider P — {Iiq, h\, /iji 5i 
the zero-dimensional ideal in W^x, u, A], generated by the polynomials defining the 



1 9m ! 9i 



System ( Sys — i I . Note that the ideal J* to be zero dimensional is equivalent to 



suppose that there is a finite number of solutions of the set of polynomial equations 
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defining the system. This set of solutions is denoted by V{J'^) and it is usually 
called the variety of J^ The goal is then to compute the points in V(J*). 

In order to do this we resort to the moment approach (see Lasserre et al. |18j . 
Lasserre [T5]). It is clear that among all the finite probability measures defined on 
the support of our feasible set T^(J*), those assigning positive probability to all the 
solutions of V{J^) give maximal range to the moment matrix M^r* (a;, it, A), and for 
sufficiently large N* rank{M]y.{x,u, A)) = |y(J')|. Moreover, since V{J^) is finite 
by Theorem |6j there is an index N* such that all the SDP relaxations SDP^ — N 



i = 1, . . . ,m, are exact and their solutions are moment sequences of probability 
measures with support on V{J''). Hence, the idea is to translate that condition, i.e. 
that we look for measures which give positive probability to all points in V{J^), 



on the moment variables of our SDP relaxations (SDP^ — N* I, i = 1, . . . ,m. In 
this way we will get the points in V{J'^) as solutions of the relaxations. This 
underlying approach was theoretically justified by Lasserre et al. |TF and Lasserre 
[161 Theorem 6.6] and can be effectively implemented by using the moment matrix 
algorithm as described in (16i Algorithm 6.1]. The reader is referred to ^181 ITS] for 
further details. 

Let us consider the quotient space M.[x,u, X]/ whose elements are the cosets 
[/] = {/ + 9 ■ 9 € J} for any / £ R[x,u,X\. Since J* is zero-dimensional, 
M[a:,u, A]/J' is a finite dimensional K-vector space with the usual addition and 
scalar product. Furthermore, M.[x, u, A]/ J' is an algebra with multiplication [/] [g] = 
[fg]. If Bji is a basis of R[x, u, A]/ J*, then the multiplication matrix M/i associated 
with the multiplication operator nih : M[a;,u, A]/ J* — > M[x,m, A]/ J*, mh{[f]) = [fh] 
for h e A], plays and important role for obtaining the points in V{J^). 

Let d = maxjmaxj^i „ Cj,max ^ ^ i...m Vl} the maximum half degree of the 

i — l...n 



polynomials defining System (Sys — i|, i = l,...,rn. Finally, let us denote by 
r (sDP' - N*\ the feasible region of Problem ^SDP' - N*\ , i = l,...,m. 

The above discussion allows us to apply to the variety V[J^), theorems 6.2 and 
6.5 in [TB| and Proposition 4.6 in [TH] which leads us to state our final result. 

Theorem 8. For N* large enough, there exists d < t < N* such that: 

rankMt(a:,M,A) = rankMt„d(a;, w. A) = \V{r)\, V(a;,u,A) e R jSDP' - N*l Vi = 1, , 



Moreover, one can obtain the coordinates of all (x, u, A) e V{J^), Vi = 1, . . . , n, as 
the eigenvalues of multiplication matrices. 



The following example, that appears in [5], illustrates the methodology proposed 
in this paper. 



Example 9. Consider problem (MOLP) with the data 



C 



A 



( 2 


1\ 






1 


1 




3 


1 


2 


, b = 


4 


-1 







-5 


^0 


-V 







Observe that the last two constraints refer to the upper bound constraints < 5 
and X2 < 5, so they are not considered as rows of the matrix A but as the sets of 
upper bounds in the polynomial constraints pj (x) in Theorem 

We use Gloptipoly 3 14 to formulate the semidefinite problems of Theorem^ 
and SEDUMI 1 . 3 ^27j as the SDP solver. Since in this problem m = 3, according 
to Theorem three semidefinite problems must be solved for each relaxation order 
N* . The original primal region has integer extreme points. Thus, M = 1 and 
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we use their original hounds for the range of values in problems (SDP^ ~ N* 
Therefore, we must only adapt the range of A. For i = 1,2,3 we get that for 
Ml = M2 = -M3 = 6 and N* = 4, the rank condition of Theorem^ is satisfied, i.e. 
rank M4(g;, u, A) = rankMi(x,/i, A) = 2, and we extract the following solutions of 
^SDP' - N*] (2 for i^l,2fori^2 and 2 for i = 3 ): 





i = 1 


i = 2 


i = 3 




X u X 


X u X 


X u X 


Sol. #1 
Sol. #2 


(1,2) (2,0,0) (4,2) 
(0,4) (2,0,0) (4,2) 


(1,2) (0,3,0) (3,3) 
(2,1) (0,3,0) (3,3) 


(2,1) (0,0,3) (3,3) 
(4,0) (0,0,2) (2,4) 



Note that the number of moments involved in the SDP problems that had to be solved 
was 1716 (923 after preprocessing /moments substitution). The moment matrix 
Mat. (x. It, A) has size 210 x 210, there are 6 localizing matrices with size 84 x 84 
and 4208 linear inequalities. 

Thus, projecting the set of extracted solutions onto the x- coordinates, we get the 
set of Pareto-optimal extreme solutions of the problem, Xe — {(4,0), (1,2), (2, 1), (0,4)}. 
These Pareto-optimal solutions and the complete Pareto-optimal set are shown in 
Fig. [7] (black dots and black segments, respectively) . 




Figure 1 . Pareto-optimal set of Example [9] 



5. Conclusions 

Due to the similarities between standard linear programming and multiobjective 
linear programming several authors wondered whether it would exist an interior 
point method valid to find the entire set of Pareto-optimal solutions of MOLP. This 
paper answers positively this question presenting an interior point based methodol- 
ogy to find that set. Our approach is constructive and we give an explicit sequence 
of m SDP problems the solutions of which are the entire set of Pareto-optimal ex- 
treme points of MOLP. Moreover, we show how all these points can be obtained by 
applying the so called moment matrix algorithm (see [13 EH]). However, although 
our construction is explicit, we do not claim that it is computationally competitive 
with other currently available methods. The main drawback is the size of the SDP 
problems to be considered which is not polynomial in the input size of MOLP. The 
interest of our results is mainly theoretical because they show, as expected, the close 
relationship between scalar and multiobjective linear programming with regards to 
the solutions techniques. In addition, our results also show the powerfulness of 
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some techniques developed in the field of polynomial optimization ([HI [I?]) to be 
applied in apparently different areas as multiobjetive optimization. 
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